Thermoelectric properties of a chain of coupled quantum dots embedded in a nanowire 
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The thermoelectric properties of a chain of coupled quantum dots (CCQD) embedded in a 
nanowire are theoretically investigated in the Coulomb blockade regime. An extended Hubbard 
model is employed to simulate the CCQD nanowire system. The charge and heat currents are 
calculated in the framework of Keldysh Green's function technique. We obtained a closed-form 
Landauer expression for the transmission coefficient of the CCQD system. The electrical conduc- 
tance (Ge), Seebeck coefflcient (S), thermal conductance, and figure of merit (ZT) are numerically 
calculated and analyzed in the linear response regime. When thermal conductance is dominated by 
phonon carriers, the optimization of ZT is determined by the power factor (pF = S^Ge)- We find 
that off-resonant tunneling processes, asymmetrical interdot electron Coulomb interactions, weak 
interdot hopping strengths and asymmetrical tunneling rates between QDs and electrodes are not 
favored in the optimization of ZT as a result they suppress the power factor of junction system. 
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I. INTRODUCTION 

Recently, considerable studies have been devoted to 
seeking efficient thermoelectric materials with the figure 
of merit (ZT) larger than 3 because there exist potential 
applications of solid state thermal devices such as cool- 
ers and power generators. The optimization of (ZT = 
S^GeT/ k) depends on the thermoelectric response fimc- 
tions; electrical conductance (Gg), Seebeck coefficient 
(S), and thermal conductance (k). T is the equilibrium 
temperature. These thermoelectrical response functions 
are related to one another. Mechanisms leading to the 
enhancement of power factor {pF = S^Ge) would also en- 
hance the thermal conductance. Consequently, it is diffi- 
cult to find ZT above 1 in conventional bulk materials.^ 
Nanotechnology development provides a possible means 
to achieve highly efficient thermoelectric materials. Re- 
cently, it has been demonstrated that ZT's of nanostruc- 
ture composites can reach impressive values (larger than 
one).* In particular, quantum dot superlattice (QDSL) 
nanowires exhibit an interesting thermoelectric property 
in that the power factor and thermal conductance become 
unrelated.* Based on this property, one can increase the 
power factor and decrease the thermal conductance si- 
multaneously to optimize ZT. A ZT value close to 2 in 
PbSc/PbTe QDSL system was reported.^ 

The reduction of thermal conductance of QDSL was 
attributed to the increase of phonon scattering rates, re- 
sulting from phonon scattering from the nanowire surface 
states and QD interface states. ^'^ Nevertheless, the rela- 
tionship between the thermoelectric response functions 
and QD quantum confinement effect including intradot 
and interdot Coulomb interactions, electron interdot hop- 
ping, and QD size fluctuation remains unclear even 
though thermoelectric properties of QDSL have been the- 
oretically studied by solving the Boltzmann equation.^ 
In ref.[9] authors considered infinite long nanowire and 
ignored electron Coulomb interactions. In the present 
study, we investigate the thermoelectric properties of fi- 



nite number of coupled quantum dots (CQDs) embedded 
in a nanowire and connected to metallic electrodes, which 
allow carrier injection into the CQDs. Thermoelectric 
properties of a single QD and double QDs embedded in 
a matrix connected to electrodes were previously studied 
by several efforts. However, in realistic QD junc- 
tions for thermoelectric application, one needs to con- 
sider a large number of serially coupled QDs, otherwise 
it is not easy to keep the large temperature difference 
cross junction, which is crucial in the implementation of 
high efficiency thermoelectric devices.^ 

Here we consider nanoscale semiconductor QDs, in 
which the energy level separations are much larger than 
their on-site Coulomb interactions and thermal energies. 
Thus, only one energy level for each quantum dot needs 
to be considered. An A^-level Anderson model is em- 
ployed to simulate the system as shown Fig. 1. Because 
we take into account a thick insulator matrix between two 
metallic leads (for instance Si02 or SiN), direct tunnel- 
ing currents between two electrodes are prohibited. In 
general, germanium or silicon QDs can be well growth 
in Si02 (SiN) matrix. The system illustrated in Fig. 
1 may find promising applications as thermoelectric de- 
vices. When N is infinite, it can be regarded as model 
system to clarify fundamental physics of one-dimensional 
strongly-correlated systems, which is one of the most 
challenging problems in condensate matter physics. For 
finite number of QDs, the transport properties of such 
QD molecular is complicate due to nonneglecting elec- 
tron Coulomb interactions. As a consequence, it is very 
tedious to do the optimization of ZT. Our analytical form 
of charge and heat formula provides an efficiency means 
to deal such difficulties. The derived transport formula 
still valid in the nonequilibrium regime. In this study, 
we find that the optimization of ZT value favors the fol- 
lowing conditions: (1) QDs with low energy level fluctu- 
ations, (2) QD energy levels lie above the Fermi level 
of electrodes, (3) T < tc ^ Uo, where tc, Uq, and F 
arc electron interdot hopping strength, on-site electron 
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Coulomb interaction, and tunneling rate, respectively, 
and (4) = Tji with F/, + Fjj kept constant, where 
rL(r_R) is the left (right) tunneling rate. 

II. FORMALISM 

The Hamiltonian of N coupled QDs connected to 
metallic electrodes can be described by the combina- 
tion of extended Hubbard model and Anderson model 
H = Hq + Hqjj : 

Ho = "^ekal^^ak^a +^ekbl^^bk^a (1) 

k.(T k,a 

,Nd]y^^bk^a + C.C 

where the first two terms describe the free electron gas of 
left and right metallic electrodes. a\ ^ {b\ ^) creates an 
electron of momentum k and spin a with energy ej, in the 
left (right) metallic electrode. Vk,i {i = l,N) describes 
the coupling between the metallic electrodes and the first 
(A^-th) QD. dl ^ {di,a) creates (destroys) an electron in 
the ^th dot. 

Hqd = ''^Eene^cT + ^'^Uene,ane^^ (2) 

where Ei is spin-independent QD energy level, and 
n-i.a = d''i^de,a- Notations Ui and Ui,j describe the in- 
tradot and interdot Coulomb interactions, respectively. 
tij describes the electron interdot hopping. Note that 
the nearest neighbor interdot hopping and interdot 
Coulomb interaction arc taken into account in Eq. (2). 

Using the Keldysh-Green's function technique, the 
charge current leaving the left and right electrodes in the 
steady state can be expressed by 

Jl = ^ j d^^L [Vl ie)ImGl, (e) - *G<, (e)] (3) 

and 

Jk = ^ j deTR[2fR{e)ImG^^^,{t) - *G'< ,(e)], (4) 

where F^ = Fi and Tn ~ F^r denote, respectively, 
the tunneling rates of the left electrode to the first QD 
and the right electrode to the A^th QD, which are as- 
sumed to be energy and bias independent for simplicity. 
/L(fl)(e) = l/[ei^-t^LiR))/kBTL^R-^ J] denotes the Fermi 
distribution function for the left (right) electrode. The 
chemical potential difference is given by fiL~ fJ-n. = eAV. 
Ti^(R) denotes the equilibrium temperature of the left 
(right) electrode, e and h denote the electron charge 



and Planck's constant, respectively. Obviously, tun- 
neling currents are determined by the on-site retarded 
Green's function (G'^^J ^{e)) and lesser Green's func- 
tion (G^^j^(e)). It is not trivial to solve N coupled 
QDs with electron Coulomb interactions (Ue and Uij) 
and electron interdot hopping {tij = tc). Conventional 
Hartree-Fork mean field theory can not resolve the de- 
tailed quantum pathes. Our previous approach beyond 
the mean-field approximation can obtain all Green func- 
tions in the closed form solutions for N=2 in the limit 
of tc/U <C 1.^^'^"^ Based on the same approach, we have 
demonstrated that the charge and heat currents of CC- 
QDs can be expressed as^'^ 

J = ^ 1 deTmL{e)-fnie)], (5) 

Q ^ IJ deTie)ie-EF-eAV)[h{e)-fRie)],ie) 

where T(e) = (71, Ar(e)-l-7Ar4(e))/2 is the transmission co- 
efficient. Te.j{e) denotes the transmission function, which 
can be expressed in terms of the on-site retarded Green's 
functions, even though Tejie) should be calculated by 
the on-site retarded and lesser Green's functions. The 
transmission function in the weak interdot hopping limit 
{tij <C U) has the following form, 

2xi^" Fi(e)F™^(e) 

where Im means taking the imaginary part of the function 
that follows, and 

Gl„M = Pi,™/(Mi - ni,„ - j:t,n), (8) 

where fj.i ^ e — Ei + iFi/2, Ili,m denotes the sum of 
Coulomb energies arising from other electron present 
in the first QD and its neighborhood, and F™^(e) = 
— 2/mE™jy(e), where S™7V denotes the self energy result- 
ing from electron hopping from QD 1 to QD A^ through 
channel m. For electron with spin a tunneling from the 
left electrode into level 1 of A^ serially coupled QD and 
exit to the right electrode, we have 2 x 4^~^ quantum 
pathes (or channels), since level 1 can be either empty or 
singly occupied (with spin —a) and all other levels can 
be empty, singly occupied (with spin up or down), and 
doubly occupied. For the N — 2 case, the explicit expres- 
sions of probability weights and self energies have been 
worked out in our previous paper. -'^'^ Here, we generalize 
the expressions to a CCQD system with arbitrary num- 
ber of QDs. We found that for the A^-QD system, the 
probability factors pi_m are determined by the following 
relation 

(di + biij)Y\_{ae + b(a + bta + C() ^ ^ Pl^rn, (9) 
i=2 m=l 
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where at = \-Ni^s, hi^s = Nt,a, ai = \~ Ni,,^- Ni,,„+ci, 
be,g- = Ni^g. - Cf, bt^a- = Ni^a- ~ Q. ai, bi^^, bf^^r, 
and Cf describe the probability factor for the £-th QD 
with no electron, one electron of spin a, one electron 
of spin (T, and two electrons, respectively. Nt^^r, Ni^^, 
and ct = {nt^gTii^g) denote the thermally averaged one- 
particle occupation numbers for spin a and a and two- 
particle correlation functions. We note that the sum of 
Eq. (9) equals one, which indicates probability conserva- 
tion. S™^ is given by a continued fraction 
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112, m — 



(10) 



where Ilim denotes the sum of Coulomb energies due to 
interaction of an electron entering the ^-th QD with the 
other electrons present in the CCQD in configuration to. 
Note that in the uncorrelated limit {Ui = and Uij = 
0), the above expressions reduce to the exact solution to 
the quantum transport through a linear chain within the 
nearest- neighbor tight-binding model. Therefore, our 
results are still able to illustrate the transport behavior 
in the other extreme, U/tij <C 1, since the terms we 
ignored are related to correlation of electrons in different 
levels, and it is valid when the Coulomb U terms become 
small. However, this study is restricted in the condition 
of tij/U ^ 1 throughout article. 

The average occupation numbers N£ ,j = Ni^g and 
C£ are determined by solving the on-site lesser Green's 
functions"'^^'^'^ iGf^{e), which take the following form 

■^< f\ ri(£)/L +r7^jv(£)/fl 



Thus, we have 

Ni = - 

and 



ri(e) + rrjv(e) 



&^ri(e)/L + r7;^(6)/fl 
V? r,(.)-,r™,(e) 



(12) 



Cl = - / — 2^ , r^m /,^ ^"iGi ,„(e). 



TT 



^.+1 ri(e) + r-^(6) 



(13) 

where 2_/m=4«-i+i denotes a sum over configurations ob- 
tained by the product 

JV 2x4""^ 

NisWiai + bt^^ + bt^g + ci)^ ^ Pi,™. (14) 

t=2 m=4"-i-|-l 

To calculate the tunneling current from of Eq. (3), we 
also need the retarded Green's function ,^(e), which is 
given by GJ_^(e) = GI,™(e)- 

-^iv,o- = ^N,s = and cm have the same forms as 
the above equations with the indices 1 and N exchanged. 



G^^(e) is obtained from G5^ ,„(e) by reversing the roles 
of QDs 1 and N. Namely, 

Gw,m(^) = PN,m|{^J■N " ^N,m " , (15) 



Y^m 

^N,l — 

and pN,m are determined by 

Af-l 



''N.N-l 



^iV-l.JV-2 
J5 

TT 2,1 

MW-2-njV-2,m---- ^^_n^ . 



(16) 



[aN + bNs)W_[ai + bi„ + bt^g + ci)^ ^ pw,™. (17) 



For QDs not in direct contact with leads (labeled by 
£ = 2,iV- 1), we have 



de 



^1 r,™W + r,Xe) 



-^"iG^_„(e), 



and 



2x4'~ 



E 



r-(e)/i + r- (e)/^ 



-_,^.^i r™(e)-frjp,(e) 



(18) 



(19) 

r?\(e) = -2/mS?\ and T1]j^{e) ^ -2/toE^jv are the 
effective tunneling rates for electrons from the ^-th QD to 
the left and right electrodes, respectively. The retarded 
Green's function G^^(e) is given by 



Pi.' 



(20) 



where the probability factors, pt^m are determined by Eq. 
(9) with the indies 1 and I interchanged. The self energies 
I]™;^ and are given by 



'.£-1 



yU£_l - n 



(21) 



£-l,m 



Mf-2— IIjx 



"72 — 

2,1 



and 



(22) 



7i 

JV-1,JV 



The explicit expressions of Pi^m and 11^ „i for all levels 
and all configurations of the N=3 case are given in the 
Appendix. As an example, for a five-dot CCQD with 
configuration described by pi,m = 01^20-03^4, sCs, we have 

'^l,m('^) = 
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(1 - Ni^g)b2aa3h,sCr, 



fJ-i - Ui:2 - 



M2-t^l,2- 



"72 — 

2,3 





3,4 



M4 -t/4 -2(74 5 - 



4,5 



(23) 

In the linear response regime, Eqs. (5) and (6) can be 
rewritten as 



■J - '-11^^ + '-ri^^ 
AV AT 



(24) 
(25) 



where there are two sources of driving force to yield the 
charge and heat currents. AT = — is the tem- 
perature difference across the junction. Thermoelectric 
response functions £ii, £12, £21, and £22 are evaluated 
by 



^12 = — ^ / deT{e){-^)Ep, (27) 



/:2i = ^ 1 rfer(6)(e - EF)i^)r, (28) 



and 



£2 



2T' 



/ rf6r(6)(6-£;^)(M!))^^. (29) 



Here T(e) and /(e) = l/[eM-^p)/''BT ^ j^j g^^.g evaluated 
under the equilibrium condition. 

If the system is in an open circuit, the electrochemical 
potential will be established in response to a temperature 
gradient; this electrochemical potential is known as the 
Seebeck voltage (Seebeck effect). The Seebeck coefficient 
(amount of voltage generated per unit temperature gra- 
dient) is defined as S = AV/AT = -Lu/iTCn). To 
judge whether the system is able to generate or extract 
heat efficiently, we need to consider the figure of merit, ^ 
which is given by 



ZT = 



{ZT)o 



f + Kph/Ke 



(30) 



Here Ge = CwjT is the electrical conductance and 
Ave = ((£22/7"^) — C.wS'^) is the electron thermal con- 
ductance. {ZT)o represents the ZT value in the absence 
of phonon thermal conductance, Uph- For simplicity, we 



assume Kph = K^phfiF. 
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3h 



is the uni- 



versal phonon thermal conductance arising from acous- 
tic phonon confinement in a nanowire,^''"^^ which was 
confirmed in the phonon wave guide. The expression of 



Kp/i = KphfiFs with Fs = 0.1 can explain well the phonon 
thermal conductance of silicon nanowire with surface 
states calculated by the first-principles method. The 
dimensionless scattering factor Fg arises from phonon 
scattering with surface impurities or surface defects of 
quantum dots.^ It is possible to reduce phonon thermal 
conductance by one order of magnitude when the QD size 
is much smaller than the phonon mean free paths. ^^'^^ 
Therefore, we adopt Fg — 0.01 as a fixed parameter and 
assume Fg is independent of the number of QDs and QD 



III. RESULTS AND DISCUSSION 

A. Three-QD junction 

In this section, we study the transport proper- 
ties of N=3 case, which were cxperimentally^'^'^"' and 
theoretically^^'^^ investigated in the nonlinear response 
regime to reveal the coherent and spin-dependent be- 
havior of carrier transport. In Ref. [27], the Kondo 
transport of triple QDs was investigated by using the 
slave-boson method to remove the double occupation for 
each QD. This study is restricted in the Coulomb block- 
ade regime. So far, few literatures have considered the 
spin-dependent thermoelectric properties of three cou- 
pled QDs. In Fig. 2 we plot the electrical conduc- 
tance (Ge), Seebeck coefficient (S), and electron ther- 
mal conductance (^e) as a function of gate voltage for 
various temperatures. We adopt the following physical 
parameters = Uq = GOFq, Uij = 20Fo, tij = IFq, 
Fi = F3 = F = Fo, El ^ E2 = Ef - 2OF0, and 
Ey, = Ep + 3OF0 — eVg. All energy scales are in units 
of Fq , a characteristic energy. A gate voltage is applied 
to tune the energy level of E^ such that the E3 level can 
be varied from being empty to singly occupied (see the 
inset of Fig. 2). The system with spin triple state as illus- 
trated in the inset of Fig. 2(a) is in the insulating state 
when £"3 is far above the Fermi energy. The conductance 
Ge reaches a maximum at E3 — Ep, whose magnitude 
decreases as the temperature increases. Meanwhile, Ge a 
a function of Vg has a Lorentz shape, whose FWHM is al- 
most independent of temperature as long as ksT /T > 1. 
This is referred to as the nonthermal broadening effect, 
which was also observed in the case of serially coupled 
quantum dots (SCQD).^^'^® The behavior of Seebeck co- 
efhcient (5) is illustrated in Fig. 2(b), and it is found 
that S vanishes when Ge reaches the maximum, a result 
attributed to the electron-hole symmetry. Here, holes 
are defined as missing electrons in electrodes below Ep. 
The negative sign of S indicates that electrons are ma- 
jority carries, which diffuse to the right electrode from 
the left electrode through energy levels above Ep. On 
the other hand, holes become majority carriers when S 
turns positive. Thus, the measurement of S can reveal 
the properties of resonant channels when compared with 
that of Ge. We see that S vanishes again at ksT = IFq 
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and eVg = SOFq, where E3 + Uj is lined up with Ep- 
The behavior of electron thermal conductance is sim- 
ilar to that of Ge and also shows the nonthermal broad- 
ening effect. Like the two-QD junction, the nonthermal 
broadening effect of the three- QD junction can be used 
to function as a low temperature filter. 

In Fig. 3, we show thermoelectric behaviors of a sys- 
tem with El = Ep ~ 20ro, E2 = Ep + lOFo - eVg, 
and i?3 = Ep, as illustrated in the inset of Fig. 3(a). 
Other physical parameters are the same as those of Fig. 2. 
Now the gate voltage is used to tune the level E2 such 
that the E2 level varies from being empty to singly occu- 
pied. Meanwhile, the E3 level will be depleted when E2 
is singly occupied. Although the behavior of Ge shown 
in Fig. 3(a) is very similar to that of Fig. 2(a), we note 
the FWHM of Ge in Fig. 3(a) is nearly twice as large. 
Furthermore, the nonthermal broadening effect for Kg 
disappears. This indicates that the "effective broaden- 
ing" of energy level of dot 2, which is not directly cou- 
pled to electrodes, is different from that of dots 1 and 3. 
The nonthermal broadening effect of Ge is an essential 
characteristic of resonant junction system. Once ksT is 
larger than the tunneling rates Fi = F3 = F, which is the 
broadening of energy levels of dots I and 3, the broad- 
ening of Ge depends mainly on the lifetime of the reso- 
nance, and becomes insensitive to the temperature factor 
1 / cosh^ {{e-Ep)/{2kBT)). In addition, we find more os- 
cillatory peaks of S in Fig. 3(b) as compared in Fig. 2(b). 
For example, when eVg = IOFq we have E2 = Ep, which 
matches with E3 and Ei + Uj. This resonance has very 
small probability weight and is non-observable in Gp, 
whereas it can be measured by S. The highly oscillatory 
behavior of S with respect to Vg indicates that carri- 
ers with high energies can diffuse to the right electrode 
through more resonant channels, which are far above Ep. 
This explains why the tail of Kg peak (near eVg = 20Fo 
and eVg ~ 40Fo) increases with increasing temperature. 

In Fig. 4, we plot the average occupation number Ng, 
Ge and S as functions of the gate voltage Vg (which is 
used to tune E2) for various temperatures for a thrce-QD 
junction system in the spin-blockade configuration, as il- 
lustrated in the inset of Fig. 4(b). Here, Ei ~ — IOFq, 
E2 = £;F + 10Fo-eV"g, Ej, ^ Ep-mVo, and C/^ j = IOFq. 
Other physical parameters are the same as those of Fig. 3. 
This configuration was considered in Ref. [26] within 
the framework of Master equation technique for study- 
ing spin-blockade behavior of three coupled QD in the 
nonlinear response regime. A spin-blockade can occur, 
because the second electron appears in dot 3 must satisfy 
the Pauli exclusion principle. As seen in Fig. 4(a), the E2 
level is tuned from being empty to singly occupied, while 
the E^ level remains singly occupied even in presence of 
the interdot Coulomb interactions. We observe a small 
bump of Ge near eVg = 20Fo. Although both Ei and E2 
become resonant at eVg = 20Fo, electrons in this reso- 
nance state can not tunnel to the right electrode through 
E3, due to the presence of i7j arising in dot 3. The max- 
imum Ge occuring at eVg = 30Fo is resulting from the 



resonance: Ei + Uj = E2 + U1 + U,j = E3 + U3 (described 
by pi,io), which is the Pauli spin blockade process for 
the case of three-QD junction. The maximum Ge in this 
spin singlet state is small than the maximum of Ge in 
spin triple state of Fig. 2. From the results of Figs. (2)- 
(4). we find that the maximum Seebeck coefficients are 
smaller than one. This is not preferred for the purpose 
of enhancing ZT. Such a result also implies that the fiuc- 
tuation of QD energy levels in CCQD will suppress ZT. 

In Fig. 5, we plot the occupation numbers Ni, elec- 
trical conductance Ge, Seebeck coefficient S, and elec- 
tron thermal conductance functions of gate volt- 
age at ksT = IFq. We consider identical energy levels 
with Ee ^ Eo ^ Ep + 30Fo - eVg and ^ = 6F0. All 
QD levels are tuned by the gate voltage from far above 
Ep to far below Ep. Other physical parameters are the 
same as those of Fig. 2. Because tej ~ te > ksT and 
it,j = > F, these thermal response functions (Ge, S, 
and Ke) display structures yielded by the electron hop- 
ping effect. The three peaks labeled by Vgi, Vg2, and V^3 
correspond to three resonant channels at Eq ~ \/2 tc, Eq, 
and Eq + \/2 tc, which are poles of the Green's function 
G5^„j(e) for channel m = 1, in which all three QDs are 
empty. The strengths of these peaks are determined by 
their probability weights, aia2a3. Another three peaks 
labeled by Vg4, Vgs, and Vge result from the resonances 
corresponding to channel m = 28 with pi^m = -^1ct^2(tC3, 
in which dots 1 and 3 are doubly occupied and dot 2 
occupied with one electron with spin a. (see Appendix) 
A remarkable result of Fig. 5 is the larger enhancement 
of the maximum Seebeck coefficient. This enhancement 
of the maximum S is due to the degeneracy of QD lev- 
els, instead from larger tc. Note that S is sensitive to 
the fluctuation of QD energy levels, but not to tc.^^ This 
will be further demonstrated in the N = 5 case. Unlike 
two-QD junction with identical QDs, the spectra of these 
thermal response functions for the three-QD junction do 
not possess symmetric behavior as a result of the interdot 
Coulomb interactions. The central dot (dot 2) feels the 
Coulomb interactions from both dots 1 and 2, while dots 
1 and 3 can only feel the interdot Coulomb interaction 
from the central dot. Such an effect can be observed in 
the behaviors of occupation numbers as functions of eVg. 

The electrical conductance Ge, Seebeck coefficient S, 
and figure of merit ZT are plotted as functions of de- 
tuning energy A = Ei — Ep with and without interdot 
Coulomb interaction at ksT = IOFq in Fig. 6, where we 
have adopted t^j = 6F0 and F/ — F/j = Fq. Wc note Ge 
is suppressed by the interdot Coulomb interactions. Such 
an effect becomes weak, when A increases. S is almost 
independent of Ui,j and it shows a linear dependence of 
A, roughly described by —kBA/(eT). The behavior of 
ZT as a function of A can be described by the function 
ZT= aA'^ /{T^cosh'^{A/2kBT)), where a is independent 
of T and A. When A ^ ksT, Ke becomes negligible. 
Therefore, the thermal conductance is dominated by Kthj 
which is assumed to be a linear function of T. ZT is de- 
termined by the power factor of S^Ge- The maximum 
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ZT occurs near /S.max = ^ksT and it is slightly reduced 
in the presence of Ui^j. The effect of Ui^j becomes negli- 
gible when the QD energy levels are far above Ep- In the 
following discussion, we will show that electron Coulomb 
interactions are important when Ep is above E^. 

Fig. 7 shows Ge, S, and ZT as functions of gate volt- 
age at ksT = lOFo and E/, = Eq = Ep + SOFo - eVg for 
various interdot Coulomb interactions. Other physical 
parameters are the same as those of Fig. 3. Compar- 
ing with the low-temperature results given in Fig. 5, 
we see that the six peaks of Gg at fc^T = IFg now be- 
come two broad peaks. We also note that the presence of 
interdot Coulomb interactions breaks the structure sym- 
metry (taking eVg ~ SOFq as the mid pint, where Ep sits 
between Ei and Ei + Ue with equal separation) . In par- 
ticular, ZTmax,> (for QD energy levels above Ep) is only 
slightly suppressed with increasing interdot Coulomb in- 
teractions. On the other hand, ZTmax.< (for QD energy 
levels below Ep) is significantly suppressed. The differ- 
ence in interdot Coulomb interactions between central 
dot and two outer dots leads to an artificial "QD en- 
ergy level fluctuation" which suppresses ZT. Our result 
indicates that the optimization of ZT prefers having the 
QD levels above Ep, although many studies of individ- 
ual QDs and two coupled QDs indicated that electron 
Coulomb interactions can enhance ZT when QD levels 
are below Ep.^^^^^ Authors in reference[29-32] adopted 
the mean field approximation, which truncates the high 
order Green's functions arising from electron Coulomb 
interactions. 



B. Effect of number of QDs in CCQD 

Since in realistic QD junctions for thermoelectric 
application, one needs to consider a large number of 
serially coupled QDs in order to accommodate rea- 
sonable temperature gradient across the junction, it is 
important to know the effect of number of QDs {N) on 
the thermoelectric properties of a CCQD system for a 
given tc- To clarify the effect of N on ZT, we plot Ge, 
5, Ke, and ZT as functions of temperature for various 
values of N by using the following physical parameters: 
Eg = Eo ^ Ep + 30Fo, Ue = Uo = SOFq, tc = 3Fo. 
The interdot Coulomb interactions have been turned off. 
Increasing the number of QDs would change the density 
of states (DOS) of the CCQD system from atomic limit 
to the band limit. For N varying from 2 to 5, in the 
absence of electron Coulomb interactions, we have the 
following energy spectra (under the condition tc/F ^ 1) 

N = 2:€ = Eo+tc, Eo-tc] 
N = S:e = Eo + V2tc, Eq, Eq - V2tc; 

N = 4:e = Eo + ^J{3 + V5)/2 Eo + ^ {i - Vh)/2 t^, 



Eo - V^(3 - V5)/2 Eo - v/(3 + V5)/2 U 
N = 5:e^Eo + Vst^, Eo + Eq, 

Eq — tc, Eq — Vitc- 

The separations between these resonant channels become 
small with increasing QD numbers. In addition, these 
resonant channels have different broadening. For exam- 
ple, three resonant channels of the N = 3 case occur at 
e ~ Eq + V2tc , Eq , and Eq — \/2tc , and their broadening 
widths are F/4, F/2, and F/4, respectively. As expected, 
the maximum ZT decreases with increasing N. Based on 
the results of Fig. 8, the reduction of maximum ZT is 
attributed to the reduction of Ge- Increasing N, the res- 
onant channels increases whereas probability weights of 
these resonant channels decrease. This explains why Ge 
is reduced with increasing N. Note that when Coulomb 
interactions arc turned off, Ge would become insensitive 
to N for the case of small tc and large A ~ Ei — Ep. 
The larger tc, the more important the QD number ef- 
fect becomes. In the high temperature regime, we find 
that the Seebeck coefficient is not sensitive to N [see Fig. 
8(b)], while fcg reduces significantly as N increases [see 
Fig. 8(c)]. Thus, fcpft, dominates the thermal conductance 
in the large N limit. Consequently, the behavior of ZT is 
essentially determined by the power factor (pF — S^Ge), 
and the trend of ZT with respect to increasing N is simi- 
lar to that of Ge ■ We also note that the dependence of all 
thermoelectric functions on N saturates once N reaches 
5 for the weak hopping strength considered, tc = iT^. 
Thus, it is sufficient to model the thermoelectric behav- 
iors of a CCQD with large N by using N = h. How- 
ever, for larger tc, the saturation behavior would occur 
at larger N . 

In order to find the optimization condition of ZT, we 
plot in Fig. 9 the occupation number of each QD, elec- 
trical conductance, Seebeck coefficient, and ZT as func- 
tions of gate voltage for various temperatures for a 5-dot 
CCQD. Other physical parameters are the same as the 
N = 5 case shown in Fig. 8. From Fig. 9(a) we see a 
small difference in the occupation number between ex- 
terior dots (A^i,A^5) and interior dots {N2, N3, N4) even 
though the interdot Coulomb interactions are turned off. 
This phenomena also occurs in the iV = 3 and = 4 
cases (not shown here). The occupation number fluc- 
tuation was also reported in Ref. [25] for the A^ = 3 
case by solving the Master equation. For eVg = 25Fo, 
electrons prefer to occupy the outer QDs (A^i and N^). 
For eVg = 75Tq, electrons prefer to accumulate in the 
interior QDs {N2,N3,N4). The maximum electrical 
conductance is suppressed with increasing temperature, 
while the peak width becomes wider. Such a behavior 
is no different from a single QD with multiple energy 
levels. However, we note that Gc.max does not have 
a Lorentzian shape even though each resonant channel 
has a Lorentzian shape. This is mainly attributed to the 
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formation of a "miniband" with the probability weight 
(1 - iVi,<,)(l - iV2^)(l - iV3,^)(l - iV4,^)(l - ^5,,,). Note 
that once U/j — 0, there are only two kind of proba- 
bility weights (1 - Ni^a) and Ng^g for each QD.^^This 
characteristic can also be demonstrated from the expres- 
sions for N=3 given in the appendix. When the applied 
gate voltage continue to increase, the CCQD turns into 
a Mott-insulator at half-filling Ne = Ng^^ = Ni^g = 0.5, 
resulting from Coulomb band gap.'^'^ (The electrical con- 
ductance almost vanishes for each QD with one elec- 
tron). The upper "miniband" with probability weight 
Ni^aN2sNz,sNi^gN^^g arises for electrons hopping be- 
tween energy levels at Eq + Uq. The Seebeck coefficient 
shown in Fig. 9(c) goes through zero at V^i, Vg2, and Vg3, 
respectively. At these applied gate voltages, the CCQD 
has an electron-hole symmetry. Note that the Seebeck co- 
efficient vanishes at the half-filling case. We see that the 
electrical conductance and Seebeck coefficient are very 
sensitive to QD energy levels. Due to the electron-hole 
symmetry, the Ge and ZT curves are symmetric (while 
S curve is antisymmetric) with respect to the mid point 
at eVg — 45ro. Thus, ZT curve has two maxima, one 
for QD energy levels above Ep and the other for QD en- 
ergy levels below Ep. The maximum ZT at ksT = SFq 
can reach 8 for Vg near IOFq and SOFq. When Kph domi- 
nates the thermal conductance, the maximum ZT values 
are correlated to the maximum power factor. The max- 
imum ZT occurs at neither good conducting state nor 
insulating state, because the maximum Ge (good con- 
ductance)is accompanied by a poor Seebeck coefficient, 
and vice versa. Note that the results shown in Fig. 9 are 
for the case with no interdot Coulomb interactions, i.e. 
Uij = 0. Once Utj are turned on, the spectra of Ge 
becomes somewhat complicated to analyze. Meanwhile, 
it would significantly lower the maximum ZT value for 
the peak with Ep > Eq. 

Electron interdot hopping strength tc is a key parame- 
ter in determining the bandwidth of the miniband, which 
would affect the thermoelectric properties. Fig. 10 shows 
Ge, S, {ZT)q and ZT as functions of te for various de- 
tuning energies A — Eg — Ep at ksT = lOFg. Here, 
we consider the case with on-site Coulomb interaction 
U = 60Fo and zero interdot Coulomb interaction. Note 
that {ZT)o corresponds to ZT in the absence of phonon 
thermal conductance. When tc is smaller than Fq, Ge 
increases quickly with respect to tc- This behavior can 
be explained by the fact of that the transmission factor is 
proportional to t^ for the 5-dot case when tc approaches 
zero. Once te is larger than Fq, Ge becomes almost sat- 
urated. Ge values at A = lOFg are very close to those at 
A = 20Fo. Such a behavior is also seen in the dotted line 
of Fig. 9(b), in which the QD energy level is tuned by the 
gate voltage. We notice that 5" is rather insensitive to tc 
in the weak hoping limit, tc/To ^ 1 (S* « —kBA/eT). In 
the absence of phonon thermal conductance, {ZT)q di- 
verges in the weak hoping limit, because the Lorenz num- 
ber L = Ke/{GeT) approaches zero. From Fig. 10(d), we 
see that ZT increases with increasing tc, reaching an op- 



timum value at tc ~ 2Fo and then decreases for higher 
tc- The reduction of ZT at higher tc arises from the 
faster reduction of in comparison with the increase of 
Ge- To reveal how the asymmetrical coupling between 
the QDs and the electrodes infiuences ZT, we also plot 
ZT versus Vg for the detuning energy A = SOFq for dif- 
ferent tunneling rates in Fig. 10(d). (See curves marked 
by filled triangles and diamonds) The maximum ZT is 
suppressed when the ratio of F^/F/j is far away from 1, 
while keeping the same average value, F^, + F^ = 2Fo. 



IV. SUMMARY AND CONCLUSIONS 

In summary, we have theoretically investigated the 
thermoelectric effects of CCQD embedded in a nanowire 
and connected to metallic electrodes. The length of the 
CCQD nanowire is finite, and shorter than the electron 
mean free path. Thus, the electron-phonon interaction 
effect can be ignored in this study. We have derived an 
expression for the transmission coefficient T(e) in terms 
of the on-site retarded Green's functions and effective 
tunneling rates involving electron hopping and Coulomb 
interactions, which keeps the same form as the Fandauer 
formula for charge and heat currents. Closed-form ex- 
pressions (which are correct in the limit tij/U <C 1) for 
the retarded Green's functions of a CCQD system with 
arbitrary number of QDs are derived. Such an analytical 
form is quite needed for finding the best ZT values. In 
addition, our formula still works in the nonequilibrium 
regime. 

In the linear response regime, electrical conductance 
Ge, Seebeck coefficient S, electron thermal conductance 
Ke, and figure of merit ZT are calculated for finite length 
CCQDs with QD number ranging from 2 to 5. The ther- 
moelectric properties of spin-dependent configurations 
(spin triple and singlet states) are studied. The nonther- 
mal broadening behavior of Ge is maintained in the three- 
QD case in the presence of electron Coulomb interactions, 
because such a behavior is the essential feature of reso- 
nant tunneling junction system. ^^'^^ Note that the non- 
thermal broadening effect of Ge can not be illustrated by 
using the mean-field approximation. "^"""^^ When QDs are 
not identical, these thermoelectric response functions are 
weaker, resulting from symmetry breaking in the junction 
structure. Unlike the case of double QDs,^" the interdot 
Coulomb interactions lead to considerable suppression 
of the maximum ZT value in the three-QD case when 
the degenerate QD energy levels are below Ep. This is 
caused by the symmetry breaking via interdot Coulomb 
interactions in the CCQD junction. 

We find that for a given tc, increasing the QD number 
N in the CCQD would suppress the maximum ZT value, 
and the reduction quickly saturates once N reaches 5 in 
the weak hopping limit {tc/U <C 1). However, for larger 
tc, the saturation behavior would occur at larger TV. In 
addition, we illustrate that the QD energy level fluctu- 
ation suppresses the maximum ZT by using a thrce-QD 



junction. We find that the Seebeck coefficient is insen- 
sitive to tc and Ue{Ue_j) when QD energy levels are far 
above Ep, and it can be approximated by a simple lin- 
ear expression, S = —kB^/eT, where A = — Ep 
and T is the equilibrium temperature. This character- 
istic oi S = —ksA/eT was also reported in an infinite 
Hubbard chain with narrow 

bandwidth.34-36 The feature 
of S* = Ay/ AT fa -ksA/eT can be utilized to realize 
a temperature sensor. Therefore, the CCQD system is 
more promising than the single QD junction for measure- 
ment of larger temperature difference. 

Acknowledgment 
This work was supported in part by National Science 
Council, Taiwan under Contract Nos. NSC 101-2112-M- 
008-014-MY2 and NSC 101-2112-M-001-024-MY3. 
^ E-mail address: mtkuo@cc.ncu.cdu.tw 
* E-mail address: yiachang@gate.sinica.edu.tw 



Pi, 6 = ai^2,5^3,s-;ni,6 = 112,6 = U2+Uj,n3_e = Us+Uj. 



pi,7 = 0162,563,(7; III J = Ui, 112,7 = U2 + Uj, IIsj = 2U,j 



Pl,8 = aib2,aC3;Iil,8 = C//,n2,8 = C/2 + 2C/j, Hg^g = U3 + 2Uj. 



Pi, 9 = ai62,aa3;ni,9 = ?7/, 112,9 = Ul, ^3,9 = 0. 



Pi, 10 = 0162, 0-63, <t; III, 10 = Ui, 112,10 = Ui+Uj, 113,10 = U3 



Pi, 11 = aib2,ab3,a; ^1,11 = Uj, 112,11 = Uj + Uj, U311 = Uj. 



Pi, 12 = 01^2, <tC3; III, 12 = Ui, 112,12 = Ul + 2Uj, 113,12 = U3 + UJ. 



Appendix A: Green's functions of N=3 



Pi, 13 = aiC2a3; Hi, 13 = 2C//, 112,13 = U2 + Ui, 113,13 = Uj 



In this appendix, we give the detailed expression of 
7i,3 for the N=:3 case. For simplicity, we assume the 
same intcrdot hopping constant between any two QDs, 



I.e. 



tc, and we denote C/1,2 ~ C/2,1 = Ui, U2.3 



f/3,2 = Uj, and ^ig = e - Ei + ire/2; i = 1,2,3. Note 
that r2 = 0, because QD 2 is not directly coupled to the 
electrodes. Eq. (5) becomes 



7-1,3(6) 



where 



„^iri(e) + ri':3(e) '-"^^ 



G'i.m(e) 



Pi, 



Ml - III, 



^1,3 



"1.3 - „ t2_ 

^J.2 - ii2,m - ^3_n3. 



(Al) 



(A2) 



(A3) 



Pi, 14 = aiC2b3,s; 



ni,i4 = 2Ui, n2,i4 = U2 + Ui + Uj, n3,i4 = U3 + Uj. 



Pi, 15 = aiC2b3,a; 



ni,i5 = 2Ui, U2.15 ^U2 + Ui + Uj, n3,i5 = 2Uj. 



Pi, 16 = aiC2C3; 



ni,i6 - 2Ui, U2,i6 ^U2 + Ui + 2Uj, n3,i6 = c/3 + 2Uj. 



Pi, 17 = ^1,(70203; III, 17 = C/i, 112,17 = Ui, 113,17 = 0. 



Pi, 18 = ^1, ,702^3, III, 18 = f^l, 112,18 = t^7 + t^,7,Il3,i8 = U3. 



and r™3 — — 2/mS™3 denotes the effective tunneling rate 
of the electron from the first QD energy level to the right 
electrode via channel m. The probability weights (pi,m) 
and Coulomb energies (11^, m) are given by 



Pi, 19 ^1, (702^3, o-; III, 19 = ?7l, 112,19 = C^7 + t^,7, Il3,i9 Uj. 



Pi,20 = Ni^ga2C3; 



Pi,i = 010203; Hi, 1 = 112,1 = 113,1 = 0. 



ni,20 = [/l, n2,20 = Ui + 2Uj, n3,20 = U3 + Uj 



Pl,2 = 010263,(7; IIi,2 = 0,112,2 = C/j, 113,2 = U3. 



Pi, 21 = Al, 362,(703; 



Pi, 3 = 010263,^; III, 3 = 0,112,3 = 113,3 = Uj. 



ni,21 = C/l + C/7,n2,21 - C/2 + Ui,U3,21 = Uj. 



Pl,4 = O1O2C3; ni,4 = 0, 112,4 = 2C/,7, 113,4 = C/3 + Uj 



Pi, 22 = Nl,sb2.ab3,a; 



Pi, 5 = 0162,(703; III, 5 = C^7, Il2,5 = U2, 113,5 = Uj. 



ni,22 = C/l + Ui,U2,22 ^U2 + Ui + Uj, 03,22 ^ U3 + Uj. 



Pi, 23 = Ni^ffb2,ab3,a; 

ni,23 = Ui + c//, n2,23 = c/2 + c// + Uj, 0,3,23 = 2c/j. 

Pi, 24 = A^1,S-62,<tC3; 
ni,24 = Ui+Ul,Il2,24 = U2 + Ul+2U,J,n3^24 = Uz + 2Uj. 
Pi, 25 = A^l, 5-62,(703; 
ni,25 = Ui + Ui, n2,25 = 2[//, n3,25 = 0. 
Pi, 26 = A^1,S-&2,<t&3,<t; 
ni,26 = C/l + [//, n2,26 = 2C// + iJj, n3,26 = U3. 

Pi, 27 = Ni^sb2.ah.,c,] 
ni,27 = f/i + ?7/, n2,27 = 2f// + C/j, n3,27 = Uj. 
Pi, 28 = A^1,S-62,<tC3; 

ni,28 - c/i + c//, n2,28 = 2t// + 2c/j, n3,28 = c/3 + C/j. 

Pi, 29 = A^l, (70203; 
ni,29 = C/l + 2C//, n2,29 = C/2 + 2^7/, n3,29 = Uj. 

Pi, 30 = A^i, 50263,(7; 

ni,3o = f/i+2c//,n2,3o = c/2+2C//+c/,/,n3,3o = u^+Uj. 

Pi, 31 = ^^1, (70263,(7; 

ni,3i = Ui + 2Ui, n2,3i = c/2 + 2Ui + Uj, 03,31 = 2Uj. 

Pl,32 = -^''1.(70203; 
ni,32 = Ui+2Ul,Il2.32 = U2 + 2Ul + 2Uj,n3,32 = U3 + 2UJ. 

For the retarded Green's unction of QD 2. wc have 

G?„,(e) = ^^^^^^ , (A4) 

/i2 — ll2,m ~ ^2,1 ~ ^2,3 
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where E^^ = ^^^^^ and E^3 - ^j^^ft;^. The proba- 
bihty factors, P2,m and Coulomb energies H^^m; i = 1,2,3 
are given by 

P2.1 = 020103; Hi, 1 = 112,1 = n3_i = 0. 
P2,2 = 020163,(7; ni,2 = 0,112,2 = C/j, 113,2 = U3. 
P2,3 = 020163, „; III, 3 = 0,112,3 = 113,3 = Vj. 

Via = 020103; ni,4 = 0, 112,4 = 21] J, 113,4 = C/3 + Uj. 

P2,5 = 0261,^03; III, 5 = C/l,Il2,5 C//, 113,5 = 0. 
P2,6 = 0261,(763,(7; III, 6 = C^l, 112,6 = C^/ + C/j, 113,6 = V3. 
P2,7 = 0261,^63,0-; III, 7 = C^l,Il2,7 = Ui + t/j, 113,7 = Uj. 

P2,8 = 0261,503; Hi, 8 = C/i, 112,8 = C//+2f7,/, 113,8 = Us+Uj. 

P2,9 = 026i,cr03;ni,g = ?7/, 112,9 = t^/,Il3,9 = 0. 

P2,io = a26i,cr63,5;ni,io = ?7/,n2,io = Ui+Uj,U3^io = U3. 
P2.11 = 0261, 0-63, ct; III, 11 = Ui, 112,11 = Ui+Uj, 113,11 = f7j. 

P2,12 = 0261,0-03; III, 12 = C//,Il2,12 = C// + 2f7j, 113,12 = Uz + Uj. 
P2,13 = 020103; Hi, 13 = [/i + 112,13 = 2[//, 113,13 = 0. 

P2,i4 = 020163, ct; 

ni,i4 = Ui + c//,n2,i4 = 2C// + c/j,n3,i4 = c/3. 

P2,15 = 020163,^; 

ni,i5 = c/i + c//,n2,i5 = 2Ui + c/j,n3,i5 = c/j. 

P2,16 = 020103; 

ni,i6 = c/i + [//, n2,i6 = 2C/7 + 2c/j, n3,i6 = c/3 + c/j. 

P2,17 = ^2,(70103; Hi, 17 = C//, 112,17 = C/2, 113,17 = C/j. 
P2,18 = /V2, (70163,0-; Hi, 18 = C//, 112,18 = U2 + UJ, 113,18 = U^ + U j . 
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P2,19 = N2,saib3y, Ili^ig = J7/, 112,19 = U2 + Uj, 113,19 = 2U^ 

ni,20 = Ul, n2,20 = U2 + 2Uj, U3,20 = C/3 + 2Uj. 
P2,21 = ^^2,3&l,ffa3; 

ni,2i = f/i + f//, n2,2i = U2 + Ui, n3,2i = Uj. 

P2,22 = N2,abl.,jb3,g; 
ni,22 = C/l + C//,n2,22 = C/2 + C// + C/j,n3,22 = Us + Uj. 
P2,23 = N2,abl,ab3,a; 
ni,23 = + f//, n2,23 = f/2 + C// + f/j, n3,23 = 2C/j. 
P2,24 = N2^sbl,9C3; 
ni,24 = C/l+C//,n2,24 = C/2 + C// + 2C/j,n3,24 = U3 + 2Uj. 

P2,2b = A^2,5-6i,<Ta3; 

ni,25 = 2C//, n2,25 = c/2 + Ui, n3,25 = t/j. 

P2,26 = N2,abi,„b3,g; 
ni,26 = 2C//, n2,26 ^U2 + Ui + Uj, n3,26 = C/3 + Uj. 



P2,27 = A^2,ff&l,<Tfe3,<T; 
ni,27 = 2C//, n2,27 = C/2 + C// + C/,7, n3,27 = 2C/j. 
P2,28 = A^2,ff&l,<TC3; 
ni,28 = 2C//, n2,28 = U2 + Ui + 2Uj, n3,28 = C/3 + 2[/j. 
P2,29 = N2^aCia3; 
ni,29 = C/i + 2C//,n2,29 = C/2 + 2C//,n3,29 = C/j. 
P2,30 = A^2,S-Ci63,g.; 

ni,3o = c/i+2c//,n2,3o = c/2+2c//+c/j,n3,3o = f/3+c/j. 

P2,31 = A^2,S-Ci53,o-; 

ni,3i = c/i + 2Ui, n2,3i = c/2 + 2?// + Uj, n3.3i = 2;7j. 

P2,32 = ^^2,(tCiC3; 
ni,32 = t/l + 2i7/,n2,32 = [/2 + 2C// + 2C/j,n3,32 = ?73 + 2C/j. 

03 „j(e) is obtained from „j(e) by exchanging the in- 
dices 1 and 3. The average occupation numbers A^^,o-, 
Niff, and C£ are determined by solving Eqs. (18) and 
(19). 
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Figure Captions 

Fig. 1. (a) Schematics for a chain of coupled quantum 
dots (CCQD) connected to metallic electrodes, where Tl 
and Tr describe the equilibrium temperature of the left 
and right electrodes, respectively, (b) The band diagram 
illustrating a CCQD connected to metallic electrodes. 

Fig. 2. Electrical conductance (Ge), Seebeck coeffi- 



cient (S) and electrical thermal conductance (we) as func- 
tions of gate voltage (Vg), which is used to tune the 
level for various temperatures. Ei ^ E2 = Ep — 20ro, 
E:i = Ef + aOFo - eVg. Ue = eOFo, Ue,j = 20Fo, 
tij = IFo, and F^ = F^ = Fq. 

Fig. 3. Electrical conductance (Ge), Seebeck coef- 
ficient (S) and electrical thermal conductance (ke) as 
functions of gate voltage (V^), which is used to tune the 
E2 level for various temperatures. Ei = Ep — 20Fo, 
E2 = Ep + lOFo — eVg, and E3 = Ep. Other physical 
parameters are the same as those of Fig. 2. 

Fig. 4. Occupation number (N), electrical conduc- 
tance (Ge), and Seebeck coefficient {S) as functions of 
gate voltage {Vg), which is used to tune the E2 level 
for various temperatures. Ei = Ep — IOFq, E2 = 
Ep + lOFo - eVg, and E3 = Ep - 6OF0. Ue^j = IOFq. 
Other physical parameters are the same as those of Fig. 2. 

Fig. 5. Occupation number (TV), electrical conduc- 
tance (Ge), Seebeck coefficient (5), and electrical ther- 
mal conductance (^e) as functions of gate voltage (V^) 
at low temperature (fc^T = IFq) for QDs with identi- 
cal energy levels Ei = Ep + SOFq — eVg. Other physical 
parameters are the same as those of Fig. 2. 

Fig. 6. Electrical conductance (Ge), Seebeck co- 
efficient (S), and ZT as functions of detuning energy, 
A — Ei ~ Ep at temperature ksT ~ IOFq with and 
without interdot Coulomb interactions. 

Fig. 7. Electrical conductance Ge, Seebeck coefficient 
(S), and ZT as functions of gate voltage (Vg) at temper- 
ature ksT = IOFq for various interdot Coulomb interac- 
tions. Ei = Ep + 30Fo - eVg and F^ = F^ = Fq. 

Fig. 8. Ge, S, Ke, and ZT as functions of tempera- 
ture for CCQDs with the QD number varying from 
2 to 5. Symmetrical tunneling rates F^ = F/j = Fq 
areused.Ei = Eq = Ep + 30Fo, intradot Coulomb inter- 
actions Ui = U = 3OF0, and electron hopping strengths 
te,j = tc = SFq.Fl = Ffl = F = Fq 

Fig. 9. Electron occupation number (TV), electrical 
conductance (Ge), Seebeck coefficient (S), and ZT as 
functions of gate voltage (Vg) for various temperatures. 
U = 30Fo, tc = 3Fo, and F = Fq. 

Fig. 10. Ge, S, (ZT)o, and ZT as functions of ic at 
ksT = IOFq for various detuning energies. U = 6OF0 
and F = Fq. The curves marked by filled triangles and 
diamonds in Fig. 4(d) are for A = 30Fo), but with asym- 
metric tunneling rates Tp ^ F^. 
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